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GAINING TWO DERIVATIVES ON A SINGULAR EORCE IN 
THE 2D NAVIER-STOKES EQUATIONS. 

ALEXEY CHESKIDOV AND LANDON KAVLIE 


Abstract. It has long been known, for the autonomous 2D Navier-Stokes 
equations with singular forcing in that there exist unique solutions which 

are globally in L^, a gain of one derivative. These classical techniques also show 
us that the solution is almost everywhere in . On the other hand, if the 
forcing term is in L^, it is known that the solution remains in globally, 
a gain of two derivatives. In this paper, we explore classical techniques to 
show that if the force is in for a € (—1,0), then the solution gains two 
derivatives globally. These methods break down for forces in H~^. In this 
scenario, we use a Littlewood-Paley decomposition in Fourier space to show 
that a solution which exists in at some time t must then remain in for 
a small interval of time [£, t + e). 


1. Introduction 


The autonomous 2D Navier-Stokes equations on the torus are given by 


^ j ut-\- {u- V)u — uAu -h Vp = / 

^ ^ \v-u = 0. 

Due to the pioneering work of Leray ([6]), we know that the 2D Navier-Stokes 
equations admit unique solutions u{t) with u G C'([0,T],L^) D L^{[0,T], H^) for 
each T > 0 with a forcing term / C Moreover, in both two and three 

dimensions, these solutions satisfy the energy inequality. That is, for each t > 0, 
our solution u{t) satisfies 


(2) \\u{t)\\l + iy f ||u(s)||^ids < IIU0II2 + 

Jo ^ 

With a “smoother” force / € then (1) admits unique solutions u S G([0, T], iJ^)n 
L°°([0, T], for all times T > 0. For more information on these results, see [3], 

[7], [8], among others. On the other hand, for many linear parabolic equations, like 
the heat equation, we find that u £ G([0, T], as long as u(0) £ and 

/ £ H°‘. The heat equation always gains two derivatives on the force, even with a 
singular force, such as / £ H~^. 

This gap between the Navier-Stokes equations, and linear parabolic equations 
such as the heat equation was studied by Constantin and Seregin for the Navier- 
Stokes equations [5] and later for the Fokker-Plank equations [4]. Their analysis 
involves using the modulus of continuity in physical space. They find that with a 
forcing term / £ with q > 4, the solution u remains Holder continuous with 

exponent 1 — 4/g. That is, in an sense, the function gains 2 — e derivatives, 
for any e > 0 (with “smoother forces” required to reach a two derivative gain). 
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To further bridge this gap, our analysis uses the technique of Littlewood-Paley 
decompositions in Fourier space to show that u remains in locally in time, a 
gain of two derivatives on the force / € H~^ = Specifically, we prove the 

following theorem. 

Theorem 1.1. Let u he the unique solution to (1) with u(0) := uq G . Then, 
there exists T := T{uQ,f) so that u G L°°([0,to\, H^) for all 0 <to <T. 

Intervals where u G are known as intervals of regularity for the Navier-Stokes 
equations. Using Theorem 1.1, one can prove a Leray characterization for the 2D 
Navier-Stokes equations with force / G That is, [0,oo) = Uj[aj,bj) with 

u{t) G for all t G This result is well-known for the 3D Navier-Stokes 

equations with a force f G L^. A detailed discussion of these results in the three 
dimensional case can be found in [3] or [8], among others. 

To complete our study, we explore the use of classical techniques in the interme¬ 
diate spaces where / G H°‘ with a G (—1,0). In this setting, we show that classical 
techniques work to give a global gain of two derivatives. That is, we prove the 
following theorem: 

Theorem 1.2. Let f G and u(0) := Uq G for some a G (—1,0). Then, 

there exists a solution u(t) to (1) so that u G L°“([0, oo), 

The outline of our paper is as follows: In Section 2, we introduce the classical 
setting for studying the Navier-Stokes equations, including the use of the Leray 
projection. We also recap the Littlewood-Paley decomposition of a function. In 
Section 3, we use the example of the heat equation to demonstrate our Fourier 
techniques for studying the Navier-Stokes equations in a “simpler” setting. In Sec¬ 
tion 4, we use classical techniques to prove Theorem 1.2. To estimate the nonlinear 
term in the intermediate space H°‘ for a G (—1,0), we use a Littlewood-Paley de¬ 
composition and Bony’s paraproduct ([1]) as demonstrated in Lemma 4.1. Finally, 
we apply our Fourier analysis techniques to prove Theorem 1.1 in Section 5. 

Throughout this paper, we use the usual convention in which C denotes an 
arbitrary, positive constant which may change from line to line. On the other 
hand, Cq, Ci, (72, etc. are fixed constants. 


2. Preliminaries 


2.1. The Projected Navier-Stokes Equations. As this paper expands on clas¬ 
sical theory, we will use classical notation as seen in ( [3] , [8] ) when exploring classical 
techniques. In Section 5, when employing harmonic analysis techniques, we use the 
unprojected Navier-Stokes equations (1) so that the reader may more easily track 
the elements such as derivatives. To begin, we apply the Leray projection Po- onto 
divergence-free vector fields to (1). For notational simplicity, we will assume that 
Paf = /. This gives us the projected equation 

(3) Ut + B{u, u) + vAu = f. 

Here, 

B{u, u) := Pcr{u ■ V)m 
Au := Pa{—A)u. 

The Stokes operator A is an unbounded, self-adjoint, positive definite operator 
with discrete eigenvalues 0 < < ^2 < ■ • ■ (we avoid the classical use of Xp as that 
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notation has a special designation which will be introduced in the next section). 
With corresponding eigenfunctions Wj, it is known that for u = 

Au = 

3 

Thus, we may define the fractional Stokes operator as 

:= ^ 

3 

Recall that the Sobolev space i7“(T^) for a € IR is a Hilbert space with the norm 

\neZ2 / 

where u{n) is the nth Fourier coefficient of u. Note that, using the operator A, an 
equivalent defintion of the 77“ (T^) norm is given by 

For the remainder of this paper, we will omit from our norms as it is implied. 
Moreover, we will use the standard convention 

\\ u\\p ■■= ||m||lp 

Remark 2.1. It is customary to project the spaces 77“ as well when using the 
projected equations. We will not make this distinction to more easily track the 
underlying space in which we are working. The fact that the evolution occurs in a 
projected space is coincidentally implied. 


2.2. Littlewood-Paley Decomposition. In this section, we briefly describe the 
Littlewood-Paley decomposition and the Littlewood-Paley theorem. This describes 
how to relate Sobolev norms in physical space via a particular breakdown in Fourier 
space. For more information on this theory, see, for example, the book by Chemin 
[2], among others. 

Choose a nonnegative radial function y G C“(1R^) so that 


x(0 = 


ia<i 

0, ia>i. 


Let := x(A;^ ^5) — xiO- For each g > 0, we let •= ^0- For technical 

reasons, let := x('C)- 

Given a tempered distribution vector field u on and g > 1, an integer, the 
gth Littlewood-Paley projection of u is given by 


Uq{x) := Aqu{x) := ^ u{k)<j)q{k)e'^’^'^ 
fcgZ2 


where u{k) is the fcth Fourier coefficient of u. Note that, by the Littlewood-Paley 
theorem. 





1/2 


for each u G 77® and s G R. 
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3. The Heat Equation 

To aid the reader in following our techniques for the Navier-Stokes equations in 
Section 5, we consider the n-dimensional heat equation on the torus T". 

/ { ut- uAm = / 

1 w(0) = uo. 

Suppose Mo G and / G H°‘ for some fixed a G R. Then, we will prove the 

following Theorem. 

Theorem 3.1. Let u be a solution to (4). Then, u G L°“([0, oo), 

Note that, using the linearity of the heat equation, one can further show that 
u G C'([0,r],iJ“+2) for all T > 0. 

The idea of our proof is to decompose our function using Littlewood-Paley theory. 
Then, apply Duhamel’s formula to the decomposition before summing. 


Proof. Multiply the first equation in (4) by (uq)q := Aq{Aqu) and integrate in 
space. This gives us that 

\ d f 

2^11^9(^)112 + ’^>^ l \\ Uq { t )\\2 < /<? ■ Uqdx . 

Applying Cauchy-Schwartz followed by Young’s inequality to the right-hand side, 
we get that 

|lK(t)||^ + uA2|K(f)||^<^||/,||^ 

Next, we apply Duhamel’s principle to get that 


\\nq{t)r,<\\uqm\le-'^^^^ 






1 - 


Multiplying by and summing in q yields the result. 


□ 


4. The Case for a G (-1,0) 

We begin our study of the 2D Navier-Stokes equations by investigating classical 
techniques. Here, we show that classical arguments yeild a uniform two derivative 
gain for / G for all a G (—1,0). Note that when a := 0 (/ G T^), it is 
known classically that u{t) G for alH > 0, a two derivative gain. We begin by 
estimating the nonlinear term in H°‘. 

4.1. Estimating the Nonlinear Term. 

Lemma 4.1. Let u G fl for j3 G (0,1). Then, 

\\Biu,u)\\H-fi < C\\u\\Hi-ff\\u\\m. 

Note that when /3 = 1, the estimate becomes 

\\B{u,u)\\h-i < C||u||2||u||hi. 




GAINING TWO DERIVATIVES IN THE NSES. 


5 


This is the classical estimate which is obtained using H’older’s inequality and the 
Ladyzhenskaya inequality. On the other hand, when /3 = 0, we may use interpola¬ 
tion to say that 

This is the classical estimate which is obtained using H’older’s inequality followed 
by Agmon’s inequality. Thus, our estimate accurately generalizes the classical 
estimates. 


Proof. Let v G . We must estimate the integral 

/ u • Vm • vdx. 


/T2 


To do so, we will use Bony’s paraproduct. Separating each term into it’s Littlewood- 
Paley pieces, we apply the necessary cancellations to find that 


u • Vu • vdx 


/T2 


— / \Up ■'^Uq ■ Vr\dx 

|p-9|<2“'T2 

r<q+l 

+ 7 / \Up ■ V Uq ■ Vr\ dx 


|p—r|<2 
q<r+l 


+ 7 / \Up ■ VUq ■ Vr\dx 


\q-r\<2 
p<r+l 

=:I + II + III. 

To estimate /, we use Holder’s inequality followed by Bernstein’s inequality to 
find that 

I< ^ ||Mp||2||Vu,||2||'yr||oo < C'||m||h-i ^ || Up || 2 Ar || Ilr || 2 ■ 

|p—9|<2 r-<p+l 

r<p+l 

Splitting the derivative Xr and applying Cauchy-Schwarz inequality gives us that 


\r-<p-|-l / \r<p-|-l 


1/2 




= .Ia 

where e <C 1 is chosen so that jd + e <1. 
For the first sum, we see that 


=:Ib 


4 = E E 

p——lr——l 


<CEA: 


2-2/3 I 

p \ 


^P\\2 


< C||u"2 


P^-1 
2 

m-0’ 
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For the second sum, 1^, we must switch the order of summation as show below: 

OO p 

p— — l r=—1 
OO OO 

r——lp—r 

OO 

<CY,\f\\Vr\\j 

r—— 1 

<C\\vfH,. 

To estimate II, we proceed as with / using Holder’s inequality followed by Bern¬ 
stein’s inequality to give us that 

II < ll^pIbllVUglloolkrlb < C||u||ffl A,||uJ|2||Ur||2- 

\p—r\<2 q<r+l 

q<r+l 

Similarly with I, we split the derivative Xq and apply Cauchy-Schwarz to get 



^ V V 

= :IIa —IIb 

Switching the order of summation in Ha and proceeding as in I, we find that 

III < C\\u\\l,., 

Ill<C\\v\\l,. 

Finally, to estimate III, we again use Holder’s inequality followed by Bernstein’s 
inequality to get 

II < Y ll“plloo||VM,||2|kr||2 < C'lkllffi Y ^pll^plblkrlb- 

\q—r\<2 p<r+l 

p<r+l 

The rest of the estimates for III proceed exactly as in the case for II. □ 

4.2. Gaining One Derivative. In this section, as well as the following section, 
we make a priori estimates. The calculations are done on the level of Galerkin 
approximations. One can then pass to the limit to obtain the stated bounds for the 
actual solutions. 

Theorem 4.2. Let f € iF“ and m(0) := uq € for some a G (—1,0). Then, 

there exists a solution u{t) to (3) so that u € L°“([0, oo), 

Proof. Taking the inner product of (3) with and integrating in space, we 

find that 

^^||u||^c+i-I-z/||u||^c +2 < [ \B{u,u) ■ A°''^^u\dx + [ \f ■ dx 

2 dt Jj2 Jj2 

< \\B{u,u)\\h'^\\u\\h<-+2 + ||/||h..||m||^.+2. 

Using Lemma 4.1 and Young’s inequality, we find that 

( 5 ) < C\\u\\h4u\\h<^^4u\\h^^^ + \\f\Uu\\H^^2 
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(6) + Hi^ii w < +kfwi 

Dropping the 77“+^ term, and using Gronwall, we have that toi t > to > 0, 

||M(7)||ffc+i < (^l|u(to)||ffc+i + ^||/||lf.(t-to)^ exp ||M(s)||^ids^ . 

After using the embedding C this becomes 

||M(i)||ffc+i < (^llM(to)llffi + - to)^ exp l|w(s)llffic^s^ • 

By the energy inequality (2), we know that 

/ l|M(s)||^ids < \\uo\\l + -||/|U-i(t - to). 

Jto ^ 

Therefore, for 0 < e < t, we have that 


\{to €[t-s,t]-. ||n(to)||Hi >M}\<^ (i||uo||l + ^ll/ll?!-) 
Letting 



we find that 

|{to e [t-e,t] : ||w(to)||ffi > M}\ < 
Therefore, there exists to € [t — e, t] so that 

So, the above argument shows us that for t > e. 



for 

Cr 

Co-.= C\\uo\\l.^^ + —\\f\\l.. 

This gives boundedness of u in for all t > e for any fixed e > 0. To 

show the boundedness of u in for small t, we go back to equation (5). Using 

interpolation, we have that 

Thus, after using interpolation and Young’s inequality on (5), we have that 


1 d 


2 

O’ 


Dropping the term, we may use nonlinear Gronwall to say that the 

norm remains bounded for small time. 

Combining this short-term bound with the previous long-term bound gives us 
that u G L°°([0, oo), as required. 
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□ 


4.3. Gaining Two Derivatives. We combine the result of the previous section 
with an analyticity argument to show the uniform gain of two derivatives in this 
regime. That is, we prove Theorem 1.2 which we restate here for the reader’s 
convenience. 


Theorem 4.3. Let f € i?“ and m(0) := Uq G for some a € (—1,0). Then, 

there exists a solution u(t) to (3) so that u G L°°([0, oo), 


To use analyticity arguments, we first need to complexify the spaces as well 
as the Navier-Stokes equations themselves. First, the complexified space is 
given by 

:= {u = ui+ iu 2 : Ml, U 2 G iF“} 
with the inner product defined via linearity as 

(mi + iU2,Vi + iV2)H^ ■= (ui,Mi)^c + (m2,M2)^c + i{{u2,Vi) - (mi,M2)^c). 

When we write (m,m)(;, we are using the usual functional pairing, in a complex 
sense. 

We will let the time t := se*®. It is known, in this setting, that there exist unique, 
analytic solutions to the Galerkin system for complex time t in some neighborhood 
of the origin. Moreover, the restriction of these solutions to the real line agree with 
the usual Galerkin approximations. 


Proof. To begin, we multiply the complexified Navier-Stokes equations by e*®, take 
the inner product with and take the real part. This gives us 

^ ^ l | w ( se *®) llff-+i + J ^ cos ( 6 »)|| u ( se *®)||^„+2 


2 ds 


= Real [e*®((B(M,M), A“+^m)(; -f (/, A“+^m)(;)] . 


Estimating the right-hand side, we first see, as in the real case, that 


|(/,A“+1m)c| < ||/||il.||m||^.+. < -^||/|||,„ 
Next, we use Lemma 4.1 to see that 

|(R(m,m),T“+^m)c| < ||B(m, M)||iL=||M||ifo+2 

< _ — _lh/lP“+® -I- 

- (l/COs( 6 l)) 2«+3 + 


-k 


i'cos{9) 

4 


hll 


2 


J2COs(d) 

4 


M 


2 


Thus, we obtain the Riccati-type inequality 



2 


j2cos(6I)||m||^„+2 < 


'cos(0) 


\\l^ 


C 

(l/COs(0))2“+3 


M 


2a+6 

i/a + l- 


This inequality shows us that for some time ||M(t)||//<»+i < M for some hxed M > 
0 and t < T ■.= T(||Mo||ifc+i, z/,/, 0). Therefore, the solutions to the complexified 
Navier-Stokes equations extend to analytic solutions in a neighborhood D of the 
origin given by 

D-.= {t = se*® :0 < s <T,\0\ < tt/2}. 

Note that D is symmetric across the real axis, by construction. Also, note that 
within D, ||M(t)||//<»+i < M. 
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Fix a compact set K C D. By Cauchy’s formula with 7 a circle in K of radius 
r < d{K, dD), we have that 



for all t G K. Taking the norm of this equation, we get that 

< -■ 
r 

Within we find that 

(7) iy\\Au\\H‘- < ||ut||Tr<^ + \\B{u,u)\\h^ + 

< + C\\u\\f-^ llullgi, + ll/llff. 

< + ll/IU. 

where we used the Poincare inequality with the Poincare constant along with 
Lemma 4.1 in the second line; we used interpolation in the second line; and we used 
Young’s inequality in the final line. 

Moving the terms to the same side of the equation and noting that 

||ut||j:/c»+i is bounded by analyticity, ||u ||2 is bounded by the energy inequality ( 2 ), 
and ||u||j:/c,+i is bounded by Theorem 4.2, we now know that ||u||j:/c <+2 is bounded 
in the compact set K . 

Using the uniform boundedness of ||u||//c<+i obtained in Theorem 4.2, we can 
rerun this argument with the same bounds at each starting point tg S [ 0 ,oo). 
Thus, ||u ||//<»+2 is uniformly bounded in a complex neighborhood of the real axis. 
In particular, \\u{t)\\H <=‘+2 < C < 00 for each t G [0,oo). □ 

5. The Case 

In this section, we present the proof of Theorem 1.1 which we again restate here 
for the reader’s convenience. 

Theorem 5.1. Let u be the unique solution to (1) with m(0) := uq G H^. Then, 
there exists T := T{uo,f) so that u G tg], for all 0 < tg < T. 

The reason analyticity methods from the previous section fails are as follows: 
Note that when a = —1, the inequality in (7) becomes 

since interpolating between the norm between and 77“+^ fails. We are now 
unable to use Young’s inequality to split the and norms to move the ||u ||//2 
terms to the left-hand side. Thus, we must use another tactic. 

Again, note that in this section, we use the unprojected Navier-Stokes equations 
so that the use of harmonic analysis techniques are more easily followed. 

Proof. Multiply the first equation in (1) with {uq)q := Aq (Aqu) and integrate in 
space. This becomes 

^^\\Uq{t )\\2 + ’^>\\\Uq{t)\\l= [ {u ■ Vu) ■ {Uq)qdx + [ fq'Uqdx. 

z at Jj 2 Jj 2 
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Apply Cauchy-Swartz and Young’s inequality to the second term on the right-hand 
side. This gives 

^^l|M,(t)||2 + i^>^l\\uqit)\\l < J^^u-Vu- {ug)gdx + ’^\\ug{t)\\l + 

For the first term on the right-hand side, use Holder’s inequality. 


/ U-VU- {Ug)q < Cllullrll 

JT2 


n|] II Ug lip 


1 1 _ 1 

T p ‘2i 

Assume that 2 < r < oo. Then, Applying the Sobolev and Bernstein inequalities 
give us that 

1 d 


ugmi + i’^>^>gmi<c\\umr\Ht)\\ 


||m9(0IIp + -^Wfqh 






For simplicity of notation, let p G (0,1) be given by p := (p — 2)/2p. Then, using 
Young’s inequality, we find that 

dt ' 






( 8 ) 


Next, we apply Duhamel’s principle. This gives 

2 


Ug(t)||^<e-‘^V||«,(0)||^ + -A-^||/g||^ 


1 — e 


-i^xU 


+ - f e’'^>-^'>Xf-^u{s)\\%,ds 
^ Jo 


Multiply this through by and sum in q. We find that 


(9) Wumrn < WMm 


+ 


c 


E 




Remark 5.2. It is worthwhile to note that inequality ( 8 ), obtained via integrating 
the force term, can also be obtained using a non-autonomous (or time-dependent) 
force. To obtain this, we need / G ^“^([ 0 , 00 ),!?“^) with a “dominating function 
in Fourier space.” We mean that there exists a p G with 

ll/,WI|2 < ||5,||2 

for alH > 0 and g > Q for some finite integer Q > —1. 


Taking a closer look at the final integral, consider the sum 

( 10 ) 

9 

Fix a constant 7 > 0 to be determined later. Then, let Qq > 1 be chosen so that 
InAg < (Ag)'*' for all q > Qq. Then, dehne 

Q{s) := minjg > Qo : < A'^P-i} 

^(■ 5 ) := ^Q(s)- 





GAINING TWO DERIVATIVES IN THE NSES. 


11 


We estimate the integral as follows 




\m < 


' [o,t]niQ(s)<Qp 


:=/ 




For /, we have that 


/< / ||u(s 
^0 


\\%^ 


^ ^ A-2 I ds 


q<Qo q>Qo 

< I I|w(s)||^i(QoAq^ + l)ds 


< c 


fMs)\\ 

Jo 


\jids. 


For //, we see that for a fixed s, (10) can be estimated by 

^ ^ A-i <Q(s)A(s)2p + l. 

q<Q(s) q>Q(s) 


By the definition of A, 2 satisfies 

22p+l^-2p-l ^ ^ ^2 ^ 8(p+ 1/2) 

But, InA < A'*' by definition. Thus, 


^{t — s) 


In A. 




8(p+l/2) 
i'(t — s) 


for 7 < 2. 

Proceding in much the same way as we did with /, we see that 
II < f l|w(s)||ffi I A(s)2p+ A’M ds 

° \q<Q{s) 9>Q(s) / 

<C /‘|Ks)|||,,(Q(s)A(s)2p + l)ds 
Jo 

<C [ ||R(s)|||,i(lnA(s)A(s)2p + l)ds 
Jo 

<C f \\u{s)\\%^{A{sf^+^ + l)ds 
Jo 

Putting these estimates together with (9), we have that 

2 „ 


Mt)rH^<\\uorH^^^nj\iH^-2 

C /■* 

H— 


1 




1 |Ks)||^.ds. 
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Next, let p := 1/4 and 7 := 1/2. Then, {2p + 7)/(2 — 7 ) = 2/3. Note that this 
choice of 7 means that Qo = 2. This gives us that 



An application of nonlinear Gronwall leads to the desired result. 


□ 
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